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Abstract 

We construct approximate formulas for the O(af) QCD contributions to vector, 
axial-vector, scalar and pseudo-scalar quark current correlators, which are valid 
for arbitrary values of momenta and masses. The derivation is based on conformal 
mapping and the Pade approximation procedure and incorporates known expansions 
in the low energy, threshold and high energy regions. We use our results to estimate 
additional terms in these expansions. 
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1 Introduction 



Correlators of heavy quark currents in different kinematical regions are of in- 
terest for a number of phenomenological applications. These two-point func- 
tions only depend on two scales, namely the square of the external four- 
momentum q 2 and the heavy quark mass m. Many of the applications focus on 
one of three distinct kinematical regions: The low energy region with q 2 « 0, 
the quark pair production threshold at q 2 = 4m 2 and the (euclidean) high 
energy region — q 2 — ► oo. 

Moments in the low energy expansion can be used for precise extractions of 
charm and bottom quark masses via sum rules pQ (for reviews see [2f3ll4"] ). 
whereas threshold and high energy expansions are directly related to produc- 
tion cross sections for ti, charmed hadrons or bottom hadrons in the respective 
energy regions. 
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The aforementioned expansions for the correlators in the three regions are 
known relatively well: In the low energy region at 0(a 2 ) the leading eight 
coefficients were computed more than ten years ago [5"|l6"fT] , and as of today as 
many as 30 moments are known [5115] . At 0(a^), however, only the first phys- 
ical moment proportional to {q 2 ) 1 for axial- vector and scalar correlators [10] 
and the first and second moment for the pseudo-scalar and vector correlators 
[T0~fllfl2|rr3"] have been available until recently. The third moment for all four 
current correlators and the fourth moment for the pseudo-scalar correlator are 
computed in Ref. [13] . 

Threshold expansions for correlators are expansions in the small heavy quark 
velocity v = yl — 4m 2 /q 2 <C 1. Currently all of the necessary machinery for 
NNLO threshold expansion is known (see for instance Ref. [15] and references 
therein). This means, that all terms of order (a™/t> n ~ 1 ) ■ {1, v, v 2 } are in prin- 
ciple known for all correlators. Explicit expansions for the axial- vector, scalar 
and vector correlators can be derived from Refs. rT5|17lll8j (see Appendix |A| . 

For high energies the leading five coefficients in the expansion are known at 
0(a 2 s ) for the case of scalar and pseudo-scalar currents [19] while for the 
vector and axial- vector currents seven terms are available [20.21] . At order 
the first two terms in the high energy expansion of the vector correlator 
have been published in Refs. [22.23]. More information is available for the 
absorptive parts of the correlators, which correspond to the logarithmic terms 
in the high energy expansions. Here the first three coefficients are known for 
the vector and axial- vector correlators [24, 25, 26 . 27128] . The leading term in 
the absorptive part has been computed recently for the vector and scalar 
correlators [23.29.30j even at order a*. 

Still, it would be desirable to have results for the correlators which are valid 
for arbitrary energies in addition to these expansions. One obvious benefit 
would be the possibility to expand such a result in a kinematical region of 
interest in order to obtain even more coefficients in the expansion. It would 
also be possible to predict values of cross sections for intermediate regions 
between threshold and high energies, where the mere expansions may not 
be very accurate anymore. Last but not least, the full energy dependence is 
essential for those QCD sum rule approaches to quark mass determination 
which use either the Borel transformation of the correlator pQ or so-called 
contour-improved perturbation theory [31|l32f 33]. 

Unfortunately, analytical results which are valid for arbitrary energies only 
exist up to 0(a s ) [51j . Still, it is possible to reconstruct the full energy de- 
pendence approximately from the known expansions at higher orders. Using 
a seminumerical method based on Pade approximations [35,36,37,38,39], the 
correlators of all four currents were reconstructed at 0(a 2 ) [5fTf4T)] . Moreover, 
it was demonstrated in [H] that in spite of the rather low amount of infor- 
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mation available at 0(az^) it is still viable to reconstruct the vector correlator 
and predict expansion coefficients with decent accuracy. Recent computations 
of additional terms in the expansions in the low and high energy domains 
[T3~|m] have confirmed these predictions and render the application of the 
approximation procedure to other correlators feasible. 

In this work we use the Pade approximation method to reconstruct the vector, 
axial- vector, scalar and pseudo-scalar heavy quark current correlators at 0(a^) 
and derive approximations to previously unknown expansion coefficients in the 
low energy, threshold and high energy regions. 

This paper is structured in the following way: In Section [2] we set up our 
conventions and explain the general approximation method. Section [3] con- 
tains details of the application of the method to the different correlators and 
the choice of physically meaningful approximants. The results for the recon- 
structed correlators and the estimates for new expansion coefficients can be 
found in Section |H In Section we summarise the work and give our conclu- 
sions. 



2 Methods 



2. 1 Polarisation functions 

It is convenient to explicitly extract the Lorentz structure of the heavy quark 
current correlators and define polarisation functions U s (q 2 ), U 5 L (q 2 )'- 

(-?V + »)nV) + »ni(? 2 ) =< / dxe^oiTjfahima) (1) 

for 8 — v, a, 
q 2 U 5 {q 2 ) =i J dxe iqx (0\T] 5 {x)j 5 {0)\0) (2) 

for 5 = s,p, 

where the currents are defined as 

3l = fairf, 3% = ^TsV', f = f = i^Hs^- (3) 

The longitudinal polarisation functions Il^(g 2 ) will not be considered any fur- 
ther in this work: U1(q 2 ) is closely related to the pseudo-scalar polarisation 
function IP(g 2 ) by a Ward identity and U v L (q 2 ) even vanishes identically. We 
do not take into account singlet contributions originating from diagrams with 
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massless cut£[] and choose the normalisation 

n 5 (o) = o. (4) 

The perturbative expansions of the polarisation functions read 

II s = n(°)- 5 + c F n«^ + n^ s f^V + (^-Y + . . . , (5) 



7T \ 7T / \ 7T 

where CV = | is the quadratic Casimir operator for the adjoint representation. 
A natural variable to describe the behaviour of U 5 is given by 

q 2 

Z= ^> (6) 
where m denotes the heavy quark mass defined in the on-shell scheme. 

For the construction of approximants we need expansions in the low en- 
ergy, threshold and euclidean high energy regions, which correspond to z — 
0, 1 and — oo respectively. Around z = the expansion reads 

o oo 

n«' 5 (z) = ^ £ <?f'V\ (7) 

The renormalisation scale \i is set to m, so that the coefficients C$ ,s are simply 
real numbers. 

Around z = 1 we have 

oo 

n w.'(*) = £ E4T(i - ^ W(i - *) • (8) 

k=k l>0 2 

The lower bound ko of the sum is 3— i for the axial- vector and scalar correlators 
and 1—i for the vector and pseudo-scalar correlators. In the naive Taylor series 

(i),5 _ jAi] ' 

k — k , 



we omit the second index, i.e. K$' s = K^' 6 

2 2 



Finally, for z — > — oo we define 



^W = EE#; W(-4.), (9) 

n=0Z>0 V ^ 

where we again use the convention D$' s = D^f . 



1 The O(a^) singlet contributions and Pade approximations to them are discussed 
in Ref. BUI. 
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2.2 Pade approximation 



The Pade approximant p n ,m(%) to a function / is denned as 



(10) 



The coefficients a i} b( are usually fixed by imposing n + m + 1 conditions of 
the form 

P&O*) (11) 

where is the j-th derivative of /. As long as all constraints have this form, 
a unique solution for the coefficients ai, bi is guaranteed [35J. 

A naive application of the Pade approximation method will, however, fail for 
the functions ns % '' s (z) because contrary to the Pade approximants (Eq. (1101) ) 
they are not meromorphic everywhere in the complex plane. There are two 
major aspects of this problem, which have to be considered: 

• The functions Il^>' s (z) diverge logarithmically for z — > — oo. There are also 
logarithmic contributions at threshold. 

• There is a branch cut along the real axis starting from z — 1. 

This behaviour can obviously not be reproduced accurately by a Pade approx- 
imation. 



The first problem related to the appearance of logarithms can be cured by 
splitting n.( l >' s (z) into two parts, 

H^ S (z)=Ui% s (z) + uSf(z), (12) 

so that H^Jg S (z) is a suitable function containing all known logarithmic con- 
tributions to Uy'' s (z). In this way the problem reduces to finding an approxi- 
mation to U^f(z). 

In the case of the vector and axial-vector polarisation functions, the functions 
U^' 5 (z) with i > 1 also show the well-known Coulomb singularity at threshold 
in addition to the logarithmic behaviour. Being proportional to (1 — z)^ l ~ 1 ^ 2 
this singularity is, however, meromorphic and can be incorporated into the 
construction of the Pade approximants. 

The second problem related to the branch cut is treated in a different way: 
We map the complex plane (including its cut) onto the unit circle in such a 
way that the branch cut is mapped onto the perimeter (see Fig. CE|). This can 
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Im(;j) 



Relzi 




Fig. 1. Conformal mapping of the complex plane onto the unit circle. The points 
z = and z = 1 correspond to u = and u = 1, respectively. z±oo goes tow = — 1. 
The branch cut starting at z = 1 is mapped onto the perimeter of the circle. 

be achieved via the conformal transformation 



Auj 



uj 



(13) 



The functions H^' s (uj) are now suitable for the Pade approximation procedure. 

In the next step we fix the coefficients hi in Eq. ([TO]) . From the expansions in 
the low energy, threshold and euclidean high energy regions we know the values 
of the II«> 5 (z) and some of their derivatives for z — 0, z — 1, and z — ► — oo. 
While the information from the first two regions can be used directly to fix 
the derivatives of II (uj) at uj = and u> = 1, respectively, the expansion 
around z — > — oo is quite different from the expansion around uj = —1. 



In order to also incorporate the information from the high energy region one 
can use the following auxiliary function 



uj 



{Auj 



in— 1 



[l+UJ 



I 271 



i (i + uj) 2 i i 



d 



j=0 



(Auy j\ \ d(i/z) 



(14) 

where n is the highest known power of 1/z in the high energy expansion. 
P n {— 1) corresponds to the coefficient of l/z n , while all other terms in the ex- 
pansion of Il^ 5 (z) around z — ► — oo together with the terms from the low en- 
ergy expansion determine the behaviour of P n (uj) around uj = 0. If the expan- 
sion around z = is available up to z m the values P n (0), P^(0), . . . , P^ +m " 1 )(0) 
together with P n (— 1) can be used for the construction of Pade approximants. 
Additional knowledge about the / leading terms in the threshold expansion 



translates into P„(T) 



KM, KM, 



pv- x) (v, 



Note that the construction Eq. ffT4l in general produces half-integer powers 
of 1/z in the high energy expansion of the reconstructed function U.^}' s (z). 
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Since only terms with integer powers of 1/z may appear in the expansion, 
we explicitely require that the terms with half-integer powers vanish. These 
additional constraints do not have the form of Eq. (TIT]) . As a result, the solu- 
tion to the system of equations which determines the coefficients of the Pade 
approximant is, in general, no longer unique. 

The Coulomb singularity proportional to (1 — z)^~ 1 ^ 2 , which appears in the 
vector and the axial-vector correlator, can easily be incorporated by defining 
P n {oS) with an additional factor (1 — u) 1 ^ 1 on the right hand side of Eq. (fill) 
so that the limit uj — > 1 becomes regular. 



3 Reconstruction of polarisation functions 



In this Section we explain the details of our calculation and list the input from 
known expansions. 



3. 1 Subtractions 



As explained in Section I2T21 the first step consists of absorbing the logarithmic 
contributions in the high energy and threshold expansions into a function 



n 



(31), S, 
log 



This function must be chosen carefully in order not to introduce 



undesired additional singularities in Il^ s (z). It is very convenient to use lower 
order analytical results as auxiliary functions [41,39], namely 



n(°).«(z) 

U w ' v (z) 



167T 2 

3 

167T 2 
1 "I 



62 



(20 4 

[J + 3z 

"5 13 M 
- + — - (1-z 
.6 Qz 

2: ^1 + 22(1-; 



4(l-z)(l + 2z) 



3z 
3 + 2z 



G{z) 



G{z) + (l-z, 



1 - 16z 

6z 



G{z) 



d \ II z) 



dz 



(15) 



(16) 



with 



Kz) 



G{z) 



=6(C 3 + 4 Lia(-u) + 2 Li 3 (n)) - 8(2 Li 2 (-«) + Li 2 («)) log 



u 



-2(2 log(l + u)+log(l 
1 \oz(u) 



2z 



u 



u 



log (u) 



1-1-1 



1-1 + 1 



(17) 
(18) 
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First, we treat the logarithmic behaviour at threshold. The relevant expansions 
of the four correlators read 

(3) _ 2.63641 -27.2677 + 0.678207^ (3)j1) 
+ (-9.47414 + 0.574190n z ) 1OR( 1 ~ _) 



+ (-17.5557 + 2.37068n, - 0.0690848nf) log(l - z) 
+ (1.31710 + 0.0312341r^ + 0.00194703nf) log(l - zf 
+ (-0.630208 + 0.0763889^ - 0.00231481nf) log(l - zf 

n (3W^ _ 2-63641 + -24.9710 + 0.678207n t + R ^ p 

1 - Z y/1 - Z ° 

+ (-9.47414 + 0.574190n,) 1 ° g } 1 Z3 

V 1 — z 

+ (-10.9576 + 2.56218^ - 0.0690848nf) log(l - z) 
+ (3.23760 + 0.00345635n, + 0.001947037^) log(l - zf 
+ (-0.630208 + 0.0763889n/ - 0.00231481nf) log(l - zf 

U i3) ' a {z) = - 0.731082 log(l - z) + K^ 3) ' a + 0{VT^) , 

n (3),s^ = _ 1 96621og(l - z) + K^ s + O(VT^z) , (19) 

where n\ is the number of light quarks. The corresponding analytic expres- 
sions and a brief outline of the derivation of the expansions can be found in 
Appendix [A] 

The construction of uf^{z) is based on the expansions of U^' v and G(z) m 
the threshold region: 

G{z) = l-±= + 0{(l-z)°), (20) 
2 Vl - z 

u (i),v ^ = _ log ^ _ z ) + const + o( v ^3^) . (21) 

Obviously, we can obtain logarithms from U^' v (z) and powers of ^= from 

Therefore, as a first attempt, we choose n^' 5 ^) to be a linear combi- 
nation of the form 

nLTW= E fci J n (1,, ^) i G( 2 ) i . (22) 

i>0,j 

This choice of U.\ g (z) will be modified later on when we consider the high 
energy region. The bounds of summation are chosen according to the powers of 
log(l — z) and } z that appear in the threshold expansions of the polarisation 
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functions. The condition that this ansatz reproduces the known logarithmic 
contributions leads to a linear system of equations for the coefficients kij, 
which determines them uniquely. 

After taking care of the threshold logarithms, we treat those appearing in the 
high energy expansions of the correlators^: 

n( 3 V(z) = - 10.0036 + 1.37572ri| - 0.0328147n 2 

+ (-0.357488 + 0.10242177, - 0.00218365n 2 ) log(-4s) 
+ (0.193107 - 0.0254675tii + 0.0004867447^) log(-4^) 2 
+ (-0.0563482 + 0.00727073^ - 0.000234540^) log(-4^) 3 

+ [ - 7.11044 + 1.01908n, - 0.031095077 2 

+ (-5.88388 + 0.75305277, - 0.014658777 2 ) log(-4z) 
+ (2.82917 - 0.25101677, + 0.00457353^) log(-4z) 2 

1 



+ 



-0.416015 + 0.034477377, - 0.000703619n 2 ) log(-4z) 3 

D (3),v + (_ 7>85787 + 0.987298/7, - 0.0260187/7 2 ) log(-4z) 

+ (0.215865 + 0.036756977, - 0.000883940n 2 ) log(-4^) 2 
+ (0.525948 - 0.0435071m + 0.000674302n 2 ) log(-4^) 3 

+ (-0.0955383 + 0.0058928177, - 0.0000879524n 2 ) log(-4^) 4 " 

U (3) ' p (z) = - 25.1130 + 3.48518r7, - 0.10285277 2 

+ (2.20686 - 0.23080877, + 0.014295777 2 ) log(-4z) 
+ (6.14249 - 0.63702477, + 0.011049977 2 ) log(-4z) 2 
+ (-1.56708 + 0.13538877, - 0.0025799477 2 ) log(-4z) 3 
+ (0.104004 - 0.0086193477, + 0.000175905r7 2 ) log(-4^) 4 

+ [ - 1.35821 + 0.177211r7, - 0.00711947r7 2 

+ (-17.9226 + 2.1094777, - 0.051567577 2 ) log(-4z) 
+ (3.44902 - 0.14487277, + 0.0011887777 2 ) log(-4z) 2 
+ (0.559407 - 0.054823777, + 0.000820889n 2 ) log(-4z) 3 

+ (-0.191077 + 0.011785677, - 0.000175905n 2 ) log(-4z) 4 



1 

I2 



In the formulas below we also give some not yet published results. Namely, the 
first two non-logarithmic terms in the 1/z expansion of the scalar, pseudo-scalar 
and axial- vector correlators come from Ref . [12] , while the logarithmic contribution 
of order 1/z 2 to the scalar and pseudo-scalar correlators come from Ref. [33] . 
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+ 



D (3), P + (_ 7 g0715 + o.685314n* - 0.0148423n 2 ) log(-4z) 

+ (-1.79679 + 0.226979n* - 0.00452222n 2 ) log(-4^) 2 
+ (1.61678 - 0.107265^ + 0.00148053n 2 ) log(-4z) 3 

+ (-0.184997 + 0.00989465^ - 0.000131929n 2 ) log(-4z) 41 ' 



iz 2 



+ o 



U (3) ' a (z) = - 9.05417 + 1.17501^ - 0.0275071n 2 

+ (-0.357488 + 0.102421n z - 0.00218365n 2 ) log(-4z) 
+ (0.193107 - 0.0254675n* + 0.000486744n 2 ) log(-4^) 2 
+ (-0.0563482 + 0.00727073^ - 0.000234540n 2 ) log(-4^) 3 

+ [15.8531 - 2.36372^ + 0.0700597n 2 

+ (-5.65395 + 0.850654^ - 0.0286846n 2 ) log(-4z) 
+ (-3.72933 + 0.420485n; - 0.00717996n 2 ) log(-4^) 2 
+ (1.15106 - 0.100911^ + 0.00187632n 2 ) log(-4z) 3 

+ (-0.104004 + 0.00861934n, - 0.000175905n 2 ) log(-4^) 4 l - 

- 1 z 

+ [D ( f' a + (6.89148 - 0.858786^ + 0.0213308n 2 ) log(-4z) 

+ (-0.537351 - 0.00471280^ + 0.000272690n 2 ) log(-4^) 2 
+ (-0.444548 + 0.0443867^ - 0.000850207n 2 ) log(-4^) 3 

+ (0.0955383 - 0.00589281^ + 0.0000879524n 2 ) log(-4z) 41 ' 
1 

U^ s (z) = - 32.1410 + 4.42783n, - 0.125454n 2 

+ (2.20686 - 0.230808n/ + 0.0142957n 2 ) hg(-Az) 
+ (6.14249 - 0.637024n, + 0.0110499n 2 ) log(-4^) 2 
+ (-1.56708 + 0.135388^ - 0.00257994n 2 ) log(-4^) 3 
+ (0.104004 - 0.00861934^ + 0.000175905n 2 ) log(-4^) 4 

+ [40.4451 - 5.58423^ + 0.173022n 2 

+ (-27.2831 + 3.37357n, - 0.0984771n 2 ) log(-4z) 
+ (-8.17872 + 0.94307072* - 0.0168387n 2 ) log(-4z) 2 
+ (4.73545 - 0.353041^ + 0.00527714n 2 ) log(-4^) 3 

+ (-0.573230 + 0.0353569n, - 0.000527714n 2 ) log(-4z) 41 1 



z 2 



+ o 



+ 



D^ s + (16.9813 - 1.84595^ + 0.0461425n 2 ) log(-4z) 
+ (-5.04041 + 0.162549^ + 0.00179429n 2 ) log(-4^) 2 
+ (-1.58168 + 0.153983ni - 0.00233074n 2 ) log(-4^) 3 

+ (0.554991 - 0.0296839n, + 0.000395786n 2 ) log(-4^) 41 ' 



z 



z 2 



+ 0[-\. (23) 
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The high energy logarithms can be generated through G(z): 

GW = z^ +0 Q). (24) 

However, the behaviour in the threshold region must not be spoiled by addi- 
tional poles induced by the threshold expansion of G(z) (see Eq. (120]) ). The 
simplest way to avoid this would be to multiply G(z) by a factor y/1 — z, which 
can however introduce non-integer powers of z in the high energy expansion. 
Following these considerations we extend the ansatz (122]) to 




where ["...] means rounding up to the next integer number. The coefficients 
d mn with m > are again fixed by demanding that the known logarithms of 
are reproduced correctly up to the highest known order. This leads to 
singular terms in the low energy expansion of nj^ ,<5 (z). We adjust the coeffi- 
cients do n in such a way that these singular terms are cancelled. Furthermore 
we can retain the property H^ s (0) = 0, which also means 11^' (0) = by 
Eq. (fT2"|) . by choosing the remaining free coefficient d 00 accordingly. 

The choice of 11^' (z) is of course not unique. A perfect reconstruction of 
the polarisation function would clearly not depend on the specific choice, so 
variations of H\Jg (z) can be used to estimate the quality of the approximation 
procedure. Following Ref. [H] , we introduce additional parameters <2j and hi 
for this purpose. For the vector and pseudo-scalar correlators we modify the 
ansatz Eq. ([25} for (z) in the following way: In the first sum, we multiply 
the summands with % = 3, j = and i = j = 1 (which roughly correspond to 
terms proportional to log 3 (l — z) and log(l — z)/y/l — z) by factors a\ + 1/z 
and a2 + 1/z, respectively. In the second sum all summands corresponding to 
the two highest powers of logarithms are multiplied by factors 1 + \f(b\z) for 
m = 3 and 1 + 1/(62-2) for m — 4. 

In the cases of the axial-vector and scalar correlator the threshold behaviour 
is quite different and only one term proportional to log(l — z) is known. Con- 
sequently the first sum in the ansatz Eq. (125]) for U log ' (z) shrinks to a single 
term. This term is multiplied by a\ + 1/z. In order to arrive at a total number 
of four parameters like in the vector and pseudo-scalar case, we multiply the 
terms corresponding to the three highest powers of logarithms in the second 
sum by factors 1 + l/(biz), 1 + l/{b 2 z) and 1 + l/(b 3 z). 

Except for the conditions eij 7^ — 1 and 6, 7^ the values of the parameters are 
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in principle arbitrary. We vary them independently with 

Oj G{-1 ± 1, -1 ± 4, -1 ± 16, -1 ± 64} , 

bt <E{±1, ±4, ±16, ±64} . (26) 

3.2 Pade approximation 

In the next step we determine the coefficients of the Pade approximants from 
the expansions in the low energy, threshold, and euclidean high energy region. 
The low energy expansions are taken from [13] ; in numerical form they read 

n (3) '"(z) =(10.6103 - 1.30278n; + 0.0282783n 2 )^ 

+ (10.4187 - 1.12407ni + 0.0223706n 2 )z 2 

+ (10.2031 - 1.01971m + 0.0194021n 2 )z 3 + 0(z A ) , 

U {3) ' p (z) =(0.812723 - 0.190853n; + 0.00721861n 2 ),2 
+ (6.33595 - 0.693155ri ; + 0.0145600njV 
+ (8.36076 - 0.803494ri; + 0.01 540 75rf)z 3 
+ (9.14377 - 0.818646ri; + 0.0149416n 2 )z 4 + 0(z 5 ) , 

n (3) ' a (z) =(4.84212 - 0.610731n/ + 0.0141353n 2 )z 

+ (2.93924 - 0.335580n ; + 0.00716845n 2 )z 2 

+ (2.06278 - 0.222971ri; + 0.00461424n 2 ),2 3 + 0(z 4 ) , 

n (3) ' s (z) =(0.123690 - 0.0679839n ; + 0.00455586nf)^ 
+ (1.78515 - 0.232769ri; + 0.00574404nf)z 2 
+ (1.92442 - 0.215014n ? + 0.00469613n 2 )z 3 + 0(z 4 ) . (27) 

The expansions around threshold and for high energies are listed in Eqs. (1191) 
and (I2"3"j) . respectively. 

3 

Following Ref. [41] . we additionally require that terms proportional to z~? 
and z~2 are absent in the high energy expansion. The resulting number of 
constraints from the different kinematic regions is listed in Table [TJ 

3.3 Discussion of approximants 

It turns out that some of the approximants have poles inside the unit cir- 
cle which translate to unphysical poles inside the complex z-plane for the 
reconstructed polarisation functions. For this reason we immediately discard 
approximants which have poles for |a>| < 1. 

There is an additional subtlety for poles which do not lie inside, but close 
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current 


low energy 


threshold 


high energy 


total 


vector 


Q 

o 


o 

z 


z+z 


y 


axial- vector 


3 





2+2 


7 


scalar 


3 





2+2 


7 


pseudo-scalar 


4 


2 


2+2 
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Table 1 

Number of constraints from the various kinematical regions for the different cur- 
rent correlators. In the high energy region constraints from both known expansion 
coefficients and absence of terms with half-integer powers of 1/z are listed 

to the unit circle. As pointed out in Ref. [U], these poles can have a huge 
numerical effect on the behaviour of the approximant on the perimeter of 
the circle. More specifically, they can lead to unphysical peaks in the cross 
sections derived from the imaginary parts of the polarisation functions above 
threshold. 

To avoid such resonances we discard approximants that have pronounced max- 
ima on the perimeter of the circle, i.e. 

>p, (28) 

where the value of p is chosen heuristically using the following two criteria: 
first, the real and imaginary parts of the polarisation functions should show 
no significant additional peaks; second, an adequate number of 0(1000) to 
(9(10000) "good" approximants should remain. In practice, we choose p = 3 
for the vector, scalar and pseudo-scalar correlators and p = 1.2 for the axial- 
vector correlator. 



4 Results 

From the Pade approximants we reconstruct the polarisation functions using 
the corresponding equations f|T2l . f|T3|) and f|T4|) . Their imaginary parts cor- 
responding to hadron production cross sections are plotted above the charm 
threshold (i.e. with n/ = 3) in Fig. [2] for all four correlators. The real parts 
are less important from a phenomenological point of view; as an example the 
vector polarisation function below and above the charm threshold is shown in 
Fig. [3j The error in the low energy region turns out to be remarkably small. 
Admittedly, this is not unexpected considering the fact that the polarisation 
functions are analytical in this region and that there is plenty of information 
from low energy coefficients. 



max 



Pn,; 



10J, 



H=i 



13 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 

v v 

Fig. 2. Imaginary part of the four loop contributions to the vector, pseudo-scalar, 
axial-vector and scalar polarisation functions above the charm threshold. The plots 
show vR®> v = vl2irIm(U^' v ), vR^> p = t>87r lm(n( 3 ) ' p ), R^< a = 127rlm(n( 3 )' a ) 
and i?( 3 )' s = 87rlm(n( 3 ) ,s ) as functions of v = y/l — 1/z. The solid black line is the 
mean from all approximants, the area covered by three standard deviations is shown 
by bands. The dashed lines show the expansions in the threshold and high energy 
regions (see Eqs. (fT9|) and ([23]) ). 

A selection of "typical" Pade approximants for all four correlators with rii = 
3, 4 and 5 massless quarks can be downloaded from 

|http : //www-ttp . particle . uni-karlsruhe . d e/Progdata/ttp09/ttp09-17/ 

The reconstructed functions can be expanded again to obtain additional low 
energy, threshold and high energy coefficients. We find that the values of the 
coefficients are strongly peaked around the mean value (for an example, see 
Fig.H]). As a consequence we give our errors in terms of standard deviations. As 
expected the error is a lot smaller for the low energy coefficients in comparison 
to the coefficients in the threshold and high energy regions. 

There are very few coefficients in the threshold and high energy regions which 
differ from the mean value by more than 50 standard deviations. We discard 
these coefficients. The resulting estimates for the expansion coefficients are 
shown in Tables [2] to O 

Comparing our predictions for the vector correlator with the previous results 
[4"T] (see Table ED, we find indications for a different sign of the threshold 
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Fig. 3. Real part of the four loop contribution to the vector polarisation function 
for the case of charm quarks. On the left hand side the region below threshold 
is shown, on the right hand side the behaviour above threshold is plotted as a 
function of v = yl — 1/ z. The dashed lines show the known expansions in the 
respective regions (see Eqs. (|19p and (|27p ). the solid lines are the mean values of 
all approximants. The shaded areas show the variation given by three standard 
deviations. In order to obtain finite values at threshold, Ren( 3 ) ,i; is plotted with an 
extra factor 1 — z below and a factor v 2 above threshold. 
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Fig. 4. Distribution of the values of C\ in on-shell scheme from different Pade 
approximants to the charm vector correlator. 



constant Kq , but good agreement for all low energy moments. Furthermore, 
the errors are decreased notably, depending on the coefficient by about one 
order of magnitude. This is mainly due to the additional information from 
Cj^' v , whereas D^' v and D^' v seem to play only a minor role. The impact 
of additional information on the quality of the predictions from the Pade 
approximants can also be seen in Table [7J where we compare results for the 
pseudo-scalar correlator with three and four low energy coefficients as input. 
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Table 2 

Expansion coefficients from the reconstructed vector and pseudo-scalar polarisation 
functions for different numbers of light quarks in the on shell scheme. c[ 3 }'^ and 

c[ 3 ]f are known exactly. The errors always apply to the last digits, i.e. 2.0(42) 
means an error of 4.2. 
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Table 3 

Expansion coefficients from the reconstructed axial-vector and scalar polarisation 
functions for different numbers of light quarks in the on shell scheme. The coefficients 
(7^3° and C^ 3 }^ are known exactly. 
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Table 4 

Coefficients from the low energy expansion of the reconstructed vector and pseudo- 
scalar polarisation functions for different numbers of light quarks. The coefficients 
are given in the MS scheme with ji = m. C^}^ and Of^f are known exactly. 





ni =3 


ni =4 


m — 5 


<?(3),a 


-2.4297 


-2.6606 


-2.7958 


C (3),a 


-3.8059 


-2.8384 


-1.9120 


G (3),a 


-2.6066 


-1.7770 


-0.9920 


C (3),a 


-1.7688(29) 


-1.1387(42) 


-0.5498(39) 


G (3),a 


-1.144(9) 


-0.692(13) 


-0.278(12) 


G (3),a 


-0.678(17) 


-0.376(24) 


-0.109(22) 


C (3),a 


-0.344(25) 


-0.166(36) 


-0.022(33) 


5 (3),a 


-0.121(34) 


-0.047(48) 


-0.004(44) 


C (3),a 


0.009(42) 


-0.004(59) 


-0.046(54) 


M3),a 
°10 


0.061(49) 


-0.024(69) 


-0.138(63) 





ni — 3 


n ( =4 


ni = 5 


C (3), s 


-5.4135 


-7.0456 


-8.1981 


g(3),s 


-12.9598 


-11.6485 


-10.3292 


q(3),s 


-6.6011 


-5.4063 


-4.2477 


q(3),s 


-3.972(14) 


-3.002(14) 


-2.073(13) 


q(3),s 


-2.665(43) 


-1.903(41) 


-1.181(39) 


q(3),s 


-1.843(81) 


-1.265(76) 


-0.723(73) 


Q(3),s 


-1.26(12) 


-0.84(12) 


-0.45(11) 


q(3),s 


-0.82(16) 


-0.54(15) 


-0.29(15) 


Q(3),s 


-0.51(20) 


-0.35(19) 


-0.22(18) 


M3),s 
U 10 


-0.30(24) 


-0.25(22) 


-0.22(22) 



Table 5 

Coefficients from the low energy expansion of the reconstructed axial-vector and 
scalar polarisation functions for different numbers of light quarks. The coefficients 
are given in the MS scheme with /x = m. The coefficients c[% a and c[% s are known 
exactly. 

5 Conclusion 



We have used the Pade approximation method to reconstruct the full en- 
ergy dependence of heavy quark correlators for vector, axial- vector, scalar and 
pseudo-scalar currents at order a^. As input we have used information from 
expansions in the low energy, threshold and high energy regions. Expanding 
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Table 6 

Comparison of the MS low energy coefficients of the charm vector correlator with 
the previous results from [41J. In the second row, the maximum error is estimated, 
whereas the error in the third row is given in terms of standard deviations. 
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Table 7 

Comparison of results from different numbers of low energy coefficients. Shown 
are the predictions for the MS low energy coefficients of the n\ = 3 pseudo-scalar 
polarisation function with three and four moments used as input. 



the reconstructed correlators, we have obtained predictions for additional co- 
efficients in these expansions. We find that these predictions are fairly accurate 
for low energy coefficients but less precise for the coefficients in the threshold 
and high energy expansions. 
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A Threshold expansion 

In this Appendix we summarise the threshold behaviour of all polarisation 
functions used in our paper. We derive the threshold behaviour of H^< s using 
Non-Relativistic QCD (NRQCD) [4~4~]|4"5] . For the present paper it is sufficient 
to obtain the results in an expansion in (1—z) up to (and including) 0((1— z)°). 
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To this end we need the matching coefficients c$ for the relation j 5 = csj^ R 
between QCD and NRQCD up to 0(a 2 ) for vector and pseudo-scalar correla- 
tors, but only up to O(a s ) for axial- vector and scalar correlators. The 0(a 2 ) 
matching coefficients are derived in Ref. jlfjflT] for the vector current (see 
also Refs. [48 49J), and in Ref. [50J for the axial- vector current. The one- loop 
matching coefficients are known since long from standard QCD computation 
(see e.g. [50J and references therein). Heavy quark current correlators reduce 
to the correlators expressed in terms of NRQCD currents, 

• J dxe^{Q\ T/(;r)/(0)|0> = c 2 i J dxe^(0\ Tj^(aOj* R (0)|0> + C 55 (q 2 ). 

(A.l) 

The second term on the right-hand side is an analytic function of q 2 , which 
corresponds to the hard heavy quark loop shrunk to a point from a diagram- 
matic point of view. This term does not contribute to R(q 2 ) oc Imil(g 2 ) due 
to its analyticity. For this reason it was never calculated in NRQCD and the 
constant Kq in the threshold expansion (jSj) remains unknown. 

For the vector case the calculation of the correlator in NRQCD was done 
analytically in Ref. [5TlfT8] . The expansion in (1 — z) can be easily obtained 
from it. For the other correlators we did not find the explicit results in the 
literature, thus we performed the calculation for the present work and present 
the result in this Appendix. 

As an illustration we show several steps for the derivation of the threshold 
expansion for the case of the vector correlator. Other correlators can be ob- 
tained in a similar way. The QCD vector current can be matched to the one 
in NRQCD by 

ft(x) = e 2 ^ (c, + ^-id ) i^a lX ](x), i = l,2,3, (A.2) 

where c\ is the matching coefficient for the vector current whose explicit form 
can be found in the references mentioned previously, and d v = 1 at the order 
of interest. The time component of the vector current vanishes in the rest 
frame of the heavy quarks with momentum q = (q , 0). Substituting the QCD 
current by the NRQCD current the right-hand side of Eq. (1A.1I) (neglecting 
C S s(q 2 )) is given by 

Cl - ^-Ej i J dxe iEx °(0\T[ X ^^](x) [^x](0)|0> , (A.3) 

where E = q — 2m = \fq 1 — 2m and integration by parts is used to relate 
the derivative to E. It is well known that the correlators in NRQCD can be 
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mapped onto Green's functions in quantum mechanics; thus one obtains 



i / dxe iExo (0\T[x^a^](x) [VV iX ](0)|0) = 2iV c %G(0, 0; E), (A.4) 



where N c = 3 is the number of colours for QCD and Sij is Kronecker's delta. 
The NNLO result for the Green's function is presented in Eq. (A. 29) of Ref. 
[T8] . One should expand the (NNLL) Green's function presented in that refer- 
ence and retain relevant terms of interest. Furthermore UV divergences should 
be subtracted according to the MS scheme conforming with the definition of 
c\ in order to arrive at the threshold expansion presented below. 



Conventionally NRQCD computations are done using the effective coupling 
(//) where the heavy quark is integrated out. We present our results using 
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(ni) 



the pole mass m and the QCD coupling constant a^ ni+1 \m) at the scale of 
the pole mass. Hence we re-express the effective coupling using the decoupling 
relation [52153] 
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with = —7/24. This induces a change in the vector and pseudo-scalar 
correlator in the coefficient of log(l — z) at order 0((1 — z 
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the vector and pseudo-scalar current correlators are given by 
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(A.6) 



Except for the term proportional to log(l —z) the result for the vector current 
agrees with Ref. |41| 3 1 . The threshold behaviour for the pseudo-scalar current 
is new. 

The scalar and axial- vector correlators are suppressed by (1 — z) compared 
to those for vector and pseudo-scalar cases. However, their leading term is 



3 The overall numerical difference is very small. For a s parametrised in terms of 
the MS heavy quark mass, we find full agreement with Ref.|41j. 
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divergent, and its 1/e pole generates a log(l — z) in NRQCD. The result reads 

n^ s (z) = - £ io g (i - Z ) + 4 3) > s + o(VT^) , 

U^ a (z) log(l - z) + ^ 3) ' a + 0(VT^~z) . (A.7) 

The origin of the log(l — z) is a triple insertion of Coulomb gluons which is of 
order O(o%) in NRQCD. The results in Eq. (I A. 71) are in agreement with the 
one in Ref. [16J. 
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